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Abstract.  The  Boltzmann  equation  for  a  gas  of  smooth,  inelastic  hard  spheres  is  introduced  and  its  homogeneous  solution 
for  an  isolated  system  is  discussed.  The  possibility  of  hydrodynamic  excitations  is  explored  using  the  linearized  Boltzmann 
equation  for  small  spatial  perturbations  of  the  homogeneous  state.  It  is  shown  that  the  spectrum  of  the  generator  for  this  linear 
dynamics  contains  five  points  at  long  wavelengths  that  are  the  origin  of  hydrodynamic  excitations.  With  this  knowledge, 
response  functions  are  introduced  that  allow  the  formal  derivation  of  linear  hydrodynamic  equations  and  associated  Green- 
Kubo  expressions  for  the  transport  coefficients  at  Navier-Stokes  order. 


1.  INTRODUCTION 

One  of  the  most  important  characteristics  of  a  real  granular  gas  is  hard  binary  inelastic  collisions.  This  feature  of 
inelasticity  alone  is  responsible  for  many  of  the  qualitative  differences  between  granular  and  normal  gases.  A  simple 
model  to  study  these  differences  is  therefore  a  system  of  N  smooth,  hard  spheres  undergoing  inelastic  collisions.  With 
the  collision  rule  specified  (see  below)  the  trajectories  in  the  N  particle  phase  space  are  well-defined  and  the  formal 
structure  of  non-equilibrium  statistical  mechanics  can  be  applied  [1,  2,  3],  In  particular,  the  BBGKY  hierarchy  for 
reduced  distribution  functions  can  be  obtained.  At  sufficiently  low  density  small  parameter  expansions  developed  for 
normal  gases  can  be  applied  for  an  approximate  closure  of  the  hierarchy  equations  [4, 5].  To  leading  order  in  the  ratio  of 
the  diameter  of  the  spheres  to  the  mean  free  path,  the  Boltzmann  equation  for  the  single  particle  reduced  distribution 
function  is  obtained1.  The  objective  here  is  to  explore  some  of  the  consequences  of  inelasticity  for  solutions  to  the 
Boltzmann  equation.  Specifically,  the  existence  of  hydrodynamic  excitations  for  small  spatial  perturbations  at  long 
wavelengths  is  demonstrated  [6,  7]  and  explored  to  Navier-Stokes  order. 

The  Boltzmann  equation  is  introduced  in  the  next  section,  parameterized  by  a  single  scalar  constant  -  the  restitution 
coefficient  a  -  differentiating  granular  gases  (a  <  1)  from  normal  gases  (a  =  I).  For  an  isolated  system  (no  external 
driving  sources)  it  is  noted  that  there  is  no  stationary  solution  due  to  the  continual  loss  of  energy  in  binary  collisions. 
However,  there  is  a  "universal"  solution  analogous  to  the  Maxwellian  whose  time  dependence  occurs  entirely  through 
a  second  moment  of  the  distribution  (the  "temperature").  This  solution  is  special  in  the  sense  that  it  appears  that 
general  homogenous  states  relax  after  a  few  mean  collisions  towards  this  solution.  This  is  similar  to  the  rapid 
velocity  relaxation  to  a  Maxwell  distribution  in  normal  gases.  A  great  deal  is  known  about  this  special  solution, 
the  "homogeneous  cooling  solution"  (HCS),  from  both  analytic  [8,  9, 10]  and  Monte  Carlo  simulation  studies  [1 1]. 

In  section  3  small  deviations  from  the  HCS  due  to  spatial  perturbations  are  described  in  terms  of  the  linearized 
Boltzmann  equation.  It  is  shown  that  the  generator  for  this  linear  dynamics  has  a  spectrum  containing  five  points  cor¬ 
responding  to  hydrodynamic  excitations  at  long  wavelengths.  The  eigenfunctions  are  identified  in  terms  of  derivatives 
of  the  HCS.  These  eigenvalues  and  eigenfunctions  show  the  existence  of  hydrodynamic  modes  for  granular  gases. 
Next,  hydrodynamic  response  functions  are  defined  from  the  formal  solution  to  the  linearized  Boltzmann  equation. 
The  hydrodynamic  spectrum  for  these  response  functions  is  obtained  to  Navier-Stokes  order  [6,  7],  and  the  transport 
coefficients  identified  in  the  form  of  Green-Kubo  expressions  [12,  13,  6], 


1  Important  questions  remain  about  the  size  of  corrections  to  the  Boltzmann  limit  due  to  velocity  correlations,  at  any  given  value  of  this 
small  parameter. 
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2.  BOLTZMANN  EQUATION  AND  HOMOGENEOUS  COOLING  SOLUTION 


Consider  a  system  of  N  smooth  hard  particles  in  a  large  volume  V.  The  particles  move  in  a  straight  line  until  a  pair  is 
at  contact.  The  velocities  of  that  pair  are  then  changed  instantaneously  according  to  a  prescribed  inelastic  scattering 
law.  In  general  this  law  can  be  quite  complex,  including  both  normal  and  tangential  friction,  the  shape  of  the  particles, 
and  a  velocity  dependent  normal  restitution  coefficient  representing  the  distortion  of  real  grains  on  collision.  After  the 
collision  the  particles  again  move  on  straight  lines  with  the  new  velocities  until  another  pair  is  at  contact.  In  this  way 
a  deterministic  trajectory  in  the  N  particle  phase  space  is  defined.  A  probability  density  in  this  phase  space  can  be 
introduced  to  represent  incomplete  knowledge  of  the  initial  state,  and  all  the  elements  of  a  non-equilibrium  statistical 
mechanics  are  then  available  for  analysis  of  this  system  [1].  In  particular,  a  Liouville  equation  for  the  probability 
density  is  obtained  and  from  it  follows  a  BBGKY  hierarchy  for  the  ^-particle  reduced  distribution  functions.  The  first 
of  these  equations  is  the  precursor  for  the  Boltzmann  equation  for  appropriate  conditions  of  low  density. 

To  simplify  the  discussion  all  subsequent  analysis  will  be  for  smooth,  inelastic  hard  spheres  of  diameter  a.  The 
Liouville  dynamics  for  the  distribution  functions  is  the  inverse  of  that  for  the  trajectories  of  the  particles  and  involves 
’’restituting”  collisions  for  the  pair  i,  j  given  by 
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where  a  is  a  unit  vector  along  the  direction  of  the  relative  position  r  J/  =  r,  —  r  •,  and  g(/-  =  v,  —  v  ■  is  the  relative 
velocity.  The  restitution  coefficient  a  takes  the  values  0  <  a  <  1,  where  a  —  I  corresponds  to  elastic  collisions.  As 
noted  above,  it  is  the  single  parameter  differentiating  granular  gases  from  normal  gases  in  this  representation.  The 
probability  density,  f(r1,yl,t),  for  the  position  and  velocity  in  the  single  particle  phase  space  obeys  the  exact  first 
BBGKY  hierarchy  equation 


dt 
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Here  /(2)(r1,r2,v1,v2,t)  is  the  joint  probability  density  in  the  two  particle  phase  space.  It  has  been  shown  elsewhere 
[4]  that  a  dimensionless  form  of  the  BBGKY  hierarchy  can  be  expressed  in  terms  of  the  single  parameter  A  =  cr/f, 
the  ratio  of  the  ’’force  range”  to  the  mean  free  path  1=  1  /(no2).  This  parameter  is  small  at  low  density,  suggesting  an 
expansion  for  a  solution  to  the  entire  hierarchy  as  a  power  series  in  A .  In  particular,  this  gives 


/'(2)(ri,ri  ±cr,V],v2,f) 
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and  the  first  BBGKY  hierarchy  equation  to  this  order  becomes  the  Boltzmann  equation 

(J?  +  v'v)/(r,v,0=7(v|/,/) 


(4) 


where  the  nonlinear  Boltzmann  collision  operator  is  identified  as 


J(y  I  /,/)  =  V2  J dy2  j da©(a  ■  g)(cr  •  g)  (a  2f(r,y\,t)f(r,v'2,t)-f(r1,yl,t)f(r2,y2,t)) 


(5) 


This  formal  expansion  of  the  hierarchy  also  gives  an  equation  for  the  corrections  to  the  Boltzmann  equation.  Therefore 
conceptual  issues  regarding  the  validity  and  context  of  the  Boltzmann  equation  for  granular  gases  can  be  studied 
quantitatively  [14,  5].  This  topic  will  not  be  explored  here. 

An  important  difference  between  normal  and  granular  gases  is  the  absence  of  energy  conservation  in  the  two  particle 
collision.  This  is  reflected  in  the  moments  of  the  collision  operator  which  can  be  calculated  directly  from  the  above 
definitions 
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The  two  zeros  on  the  right  side  of  (4)  correspond  to  conservation  of  mass  and  momentum.  The  last  term  results  from 
non-conservation  of  energy  and  §  is  the  fractional  rate  of  decrease  in  the  energy  due  to  the  inelasticity.  The  kinetic 
’’temperature”  T  and  flow  velocity  u  are  defined  in  the  usual  way 


(tT  )=/"-(  -'"<v"U)2)/<r'V',')- 


(V) 


The  "cooling  rate"  E,  [/,/]  is  a  positive  functional  of  /.  This  implies  that  the  energy  of  an  isolated  system  decreases 
monotonically.  Consequently,  there  is  no  stationary  state  such  as  the  Maxwell-Boltzmann  distribution  for  normal  gases. 
Instead,  there  is  a  homogeneous  cooling  solution  with  the  scaling  form 

/to(vTW)  =  nv03(T(t))<l>  (v/v0(T(t)))  ,  v0(T(/))  =  y/2T{t)/m.  (8) 

An  equation  for  the  temperature  follows  directly  from  this  scaling  and  the  Boltzmann  equation  without  further 
assumption 

(4  +  £(r(0))r(0  =  o.  (9) 

where 

2 
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The  temperature  dependence  of  §  can  be  made  explicit  through  the  scaling  form  of  /,  ,  so  (9)  can  be  integrated  for 
the  exact  time  dependence  of  the  temperature.  The  function  /,  is  determined  from  the  Boltzmann  equation  which 
now  takes  the  form 

\t>  Vv  {yfhcs)  =  J  (vl  fhcs’fhcs)  ■  (H) 

This  equation  must  be  solved  self-consistently  with  the  equation  for  q  to  determine  the  HCS.  Although  an  explicit 
solution  has  not  yet  been  obtained,  its  form  at  large  and  small  velocities  has  been  determined  by  good  analytic 
approximations  [8,  9].  More  generally  the  universality  of  the  HCS  as  the  attractor  for  a  wide  class  of  initial  conditions, 
its  functional  form,  and  the  scaling  property  have  been  verified  by  Monte  Carlo  simulation  [11].  In  all  of  the  following 
the  HCS  will  be  treated  as  a  known  function. 


3.  LINEAR  BOLTZMANN  EQUATION  AND  HYDRODYNAMIC  MODES 


As  noted  above,  there  is  evidence  that  small  homogeneous  perturbations  of  the  HCS  decay  on  the  order  of  the  time  for 
a  few  collisions.  A  description  of  the  response  to  small  spatial  perturbations  of  the  HCS  is  both  more  interesting  and 
more  complex,  and  is  the  object  of  the  remaining  discussion  here.  The  simplest  case  is  an  initial  spatial  perturbation 
of  the  HCS  (other  examples  might  arise  from  small  external  forces  or  boundery  condtions).  The  distribution  function 
is  written  as 

./  (r:  v./)  ./;„,(v:v(0)  [1 +  e'k  rA(k,v,f)]  (12) 

and  solutions  are  sought  to  linear  order  in  A.  Since  only  linear  terms  are  retained,  it  is  sufficient  to  consider  only  a 
single  Fourier  mode.  For  the  moment,  questions  of  stability  are  deferred.  The  Boltzmann  equation  then  takes  a  linear 
form  for  A.  Due  to  the  time  dependence  of  the  reference  state  this  equation  depends  explicitly  on  time.  It  is  convenient, 
therefore,  to  introduce  the  dimensionless  variables 


f(\,t)  =nv03(T(t))f*(\*is), 
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where  i  =  1  /no3  is  the  mean  free  path  and  .v  is  the  average  number  of  collisions  in  the  time  interval  (().?).  The 
linearized  Boltzmann  equation  then  becomes 


(ds  +  tk*  -v*  +  A)  A*  =  0 


(14) 
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where  the  linear  collision  operator  A  is  defined  by 


A x  =  LX  +  \^fhc)y,-(yfhcsX)  (15) 

In  these  dimensionless  variables  A  is  time  independent.  The  first  term  on  the  right  is  the  linearized  Boltzman  collision 
operator 

LX  =  TfTtXVf  [J(V  1  fhcsJhcs*)  +  I  fhcsX’fhcs)\ 

'oA  \T))Jhcs 

while  the  second  term  describes  the  effects  of  cooling  explicitly,  and  arises  because  the  derivative  with  respect  to  ,v  is 
taken  at  fixed  v*  . 

The  formal  solution  to  the  linear  Boltzmann  equation  is 

A*  (kW)  =  e(,'k‘-v‘+A)*A*  (k*,v*,0) .  (16) 

This  shows  that  all  linear  excitations  of  the  granular  gas  are  determined  from  the  spectrum  of  ik  *  -v*  +  A  (the  relevant 
Hilbert  space  is  defined  in  the  next  section) 

(ik*-v*  +  A)  V„(k*)  =  ln(k)Wn{k*).  (17) 


In  particular,  if  a  hydrodynamic  description  is  applicable  to  granular  gases  the  corresponding  excitations  must  occur 
in  this  spectrum.  For  a  normal  gas  with  elastic  collisions  the  hydrodynamic  spectrum  is  defined  in  terms  of  the  k  *  =  0 
spectrum.  In  that  case  there  are  five  eigenfunctions  (the  summational  invariants)  with  zero  eigenvalue,  corresponding  to 
the  conservation  of  mass,  energy,  and  momentum.  The  hydrodynamic  spectrum  for  k  *  ^  0  is  that  which  is  continuously 
connected  to  the  five  k*  =  0  zero  eigenvalues. 

A  similar  analysis  applies  for  granular  gases.  The  hydrodynamic  eigenvalues  for  k  *  =  0  can  be  identified  by 
differentiating  the  equation  for  the  HCS  with  respect  to  the  average  density,  temperature,  and  flow  velocity  (denoted 
generically  by  yd)  [6,  7] 


d 

dyr 


2^vv  (v.4J  -  J  {y\fhcs,fhcs) 


=  0. 
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In  this  way  five  eigenfunctions  and  eigenvalues  are  obtained 


A  ({)„  =  A„(O)0„,  <t>„  =  \j/n  (0)  (19) 

^(0)^  (20) 

*7  (4  ■ +  \nfhcs,  -3  -  In  fhcs,  -  Vy,  In  f^)  .  (21) 

These  are  the  hydrodynamic  modes  since  the  eigenvalues  are  the  same  as  those  of  the  exact  macroscopic  balance 
equations  for  mass,  energy,  and  momentum,  linearized  and  taken  at  k  *  =0.  This  is  therefore  a  direct  confirmation  of 
hydrodynamic  excitations  for  a  granular  gas,  independent  of  other  derivations  such  as  the  Chapman-Enskog  expansion. 
Of  course  there  may  be  other  non-hydrodynamics  solutions  for  A  (0),  but  to  simplify  the  notation  it  will  be  understood 
below  that  A(0)  refers  specifically  to  the  set  of  five  values  above.  For  small  but  finite  wavevectors,  k  '  -f  0,  the 
corresponding  eigenfunctions  and  eigenvalues  can  be  calculated  from  (17)  by  perturbation  theory  using  (19)  -  (21) 
as  the  reference  state.  The  results  obtained  in  this  way  to  Navier-Stokes  order  agree  in  detail  with  those  from  the 
Chapman-Enskog  expansion  [6,  7].  Instead,  an  alternative  approach  based  on  response  functions  will  be  explored 
here. 


4.  HYDRODYNAMIC  RESPONSE  FUNCTIONS 

To  simplify  the  notation  it  is  useful  to  introduce  a  scalar  product 

(< a,b )  =  J  dv*f;ics{v*)a+(y*)b(y*) 


(22) 
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and  a  set  of  functions  biorthogonal  to  the  k  *  =  0  eigenfunctions 

1  *2  1 
-v  +  - 

,3  2y 

Now  consider  an  initial  perturbation  of  the  HCS  in  terms  of  the  eigenfunctions  ( 

A*  (k* ,  V* , 0)  =  0m  (v*)  <5.Vm(k,0). 

The  coefficients  8ym  (k* ,  0)  can  be  determined  by  calculating  the  average  of  0„, 

(0„,A*(O))  =  5y„(k*,0). 


=  8,h 
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Since  the  0„  are  linear  combinations  of  the  summational  invariants,  the  <5  v,,(k.  0)  are  linear  combinations  of  the  inital 
deviations  for  the  density,  temperature,  and  flow  velocity  from  their  values  in  the  HCS.  Thus  the  initial  state  is  a 
hydrodynamic  perturbation.  Their  response  at  some  later  time  is  therefore 

8yn  (k“,.y)  =  G„„,  (k* ,  s)  8ym  (k* ,  0) .  (26) 

The  response  functions  Gnm(k*,  jj  are  given  by 

Gnm{ k  ,.v)  =  ^0H,0m(k  =  ^0«  (k  ,s),0m^  ,  (27) 

where  the  time  dependent  functions  are 

0m(kV)  =  e-(/k*'v*+jV)s0m,  0„(kV)  =  e-i-** V*+A+)S0„.  (28) 

The  second  equality  of  (27)  defines  the  adjoint  operator  A +  . 

Other  response  functions  can  be  defined  for  more  general  initial  conditions  and/or  for  other  variables.  These 
response  functions  defined  here  will  be  called  ’’hydrodynamic”  since  the  variables  being  considered  are  hydrodynamic 
fields,  and  the  initial  perturbation  is  such  that  the  response  is  entirely  hydrodynamic  in  the  long  wavelength  limit 

G„m( k*  =  0..v)  =  e  Xn(0]s8nm.  (29) 

The  utility  of  this  special  feature  is  demonstrated  in  the  next  section.  In  general,  for  k  *  ^  0  these  ’’hydrodynamic” 
response  functions  will  have  both  hydrodynamic  and  non-hydrodynamic  excitations 

Gnm{  kV)  =Anm(k*)e-x"^s  +  Bnm(k*,s)  (30) 

where  Bnm  (k* ,  ,sj  denotes  the  non-hydrodynamic  contribution.  It  is  expected  that  this  part  of  the  spectrum  corresponds 
to  faster  relaxation  than  that  of  the  hydrodynamic  modes,  so  that  a  hydrodynamic  description  dominates  on  a  long  time 
scale. 


4.1.  Dispersion  relations 

The  response  functions  G„m( k* ,  ,sj  describe  all  possible  excitations  in  response  to  the  initial  hydrodynamic  perturba¬ 
tion.  Generally,  it  has  the  form  of  Eq.  (30).  These  excitations  are  defined  by  the  singularities  of  the  Laplace  transform 
of  G„m(k*  ,s)  in  the  complex  plane,  e.g., 

Gnm  (k*  >  *)  =  2^  l dzeZt  +  *  (°)  +  M(k* ,  z))™l  •  (31) 

Here  I  is  the  unit  matrix,  A (0)  is  a  diagonal  matrix  whose  elements  are  Am(0),  and  A7(k' , z)  is  a  matrix  representing 
the  effects  of  finite  wavevector.  The  integral  is  defined  as  a  contour  enclosing  all  singularities  of  the  integrand.  These 
singularities  define  the  spectrum  of  the  response  functions,  i.e.  the  solutions  z(k*)  to  the  equations 

det(z(Jt*)/  +  A(0)+M(k*,z(Jt*)))  =  0.  (32) 

These  are  the  dispersion  relations  that  characterize  the  physical  nature  of  granular  gases  on  all  length  and  time  scales. 
The  specific  case  of  hydrodynamics  corresponds  to  those  solutions  which  continuously  approach  the  above  values  A  (0) 
as  k*  —y  0.  The  hydrodynamic  excitations  therefore  can  be  identified  by  determining  the  form  of  the  matrix  M( k  *,z) 
as  a  function  of  k*  and  z.  This  is  addressed  in  the  next  subsection. 
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4.2.  Green-Kubo  Relations 

A  formal  identification  of  Af(k*,z)  is  possible  using  an  analysis  paralleling  that  for  normal  fluids  [15].  In  that  case 
the  hydrodynamic  response  functions  are  constructed  from  the  local  conserved  quantities  (summational  invariants), 
and  the  conservation  laws  are  used  to  express  M  in  terms  of  response  functions  for  the  corresponding  fluxes.  The 
special  choice  of  initial  perturbations  considered  here  in  terms  of  (j)„  leads  to  similar  ’’conservation  laws” 

(ds  +  Am(0))(^m(k*,j)  +  /k*  •  ym(k*,  s)  =  0,  (33) 

with  the  ’’fluxes”  given  by 

ym(k*,s)  =  e~('k*'v*+A^ym,  Ym=v<pm  (34) 

Similarly,  the  adjoint  ’’conservation  laws”  are  found  to  be 

(«?.,  + A, „(0))^,„(k *,s)  -  t'k-y„(k*,s)  =  <5m2v(k*,s),  (35) 

with  the  conjugate  fluxes 

y„(kV)  =  e~(~,k*v*+A+)'y„,  ym  =  v0m.  (36) 

Also,  there  is  an  energy  loss  source  term  in  (35) 

v(M)  =  e-(~^+A+>v,  v  =  (1  -P)  (1  -  a2)  || dv*2f;ws(v*2,t)g*3  (37) 

where  P  is  the  projection  operator 

PX=(X,kn)$m.  (38) 

These  conservation  laws  lead  to  corresponding  equations  for  the  response  functions,  and  expose  some  important 
leading  order  dependence  of  M( k* ,  z)  on  k* .  It  is  clear  from  (29)  that  M( k* ,  z)  vanishes  at  k*  =  0.  Assuming  analyticity, 
the  exact  forms  for  the  order  k*  and  k*2  =  0  can  be  determined.  These  give  the  Euler  and  Navier-Stokes  forms  for  the 
hydrodynamics,  respectively.  The  method  is  outlined  in  the  Appendix  with  the  result 


M{k,z )  =  ikM^  +  k2M(2\z)  +  ..  (39) 

Mill  =  (find)  (4°) 

=k^/o°°dse-^  [(y^  (1  ~P)Ym)  +  8*2 £ds  v^'y1,)]  (41) 

The  terms  1 )  at  linear  order  in  k*  are  just  the  coefficients  of  the  first  order  gradients  in  the  linearized  Euler  equations. 
The  terms  of  second  order  contain  information  about  all  finite  wavevector  excitations  through  its  dependence  on  z-  As 
discussed  above,  the  special  hydrodynamic  excitations  are  those  arising  from  the  eigenvalues  of  A  (0)  at  k  *  =  0,  i.e. 
the  solutions  to 

det  (zh{k)I  +  A(0)  +  ikM™  +  k2M <2)  (zh{ 0)))  (42) 

The  contributions  to  these  solutions  from  M1'2'1  (— z,  (0))  define  the  transport  coefficients  at  Navier-Stokes  order.  From 


Eq.  (41)  there  are  terms  that  are  time  integrals  of  the  flux  -  flux  response  function.  This  is  the  generalization  of  Green- 
Kubo  expressions  for  a  granular  gas  [13,  6],  In  addition,  the  second  term  of  (41)  gives  transport  coefficients  associated 
with  the  cooling  rate  that  have  no  analogue  for  normal  gases.  The  results  given  here  agree  with  those  obtained  by  other 
methods,  such  as  the  Chapman-Enskog  solution  to  the  Boltzmann  equation  and  a  perturbation  theory  calculation  of 
the  eigenvalues  in  Eq.  (17) 
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5.  CONCLUSIONS 


The  primary  results  here  have  been  the  identification  of  hydrodynamic  excitations  in  the  spectrum  of  the  linearized 
Boltzmann  operator,  and  using  the  corresponding  eigenfunctions  to  construct  hydrodynamic  response  functions.  While 
hydrodynamics  can  be  obtained  for  granular  gases  in  other  ways  (e.g.,  the  Chapman-Enskog  expansion),  the  approach 
via  the  spectrum  is  more  direct  and  free  of  many  conceptual  stumbling  blocks  for  the  Chapman-Enskog  method. 
The  response  functions  are  special  in  the  sense  that  they  involve  the  functions  <pn  which  obey  "conservation  laws" 
relating  them  to  the  gradients  of  fluxes.  They  therefore  have  a  close  analogy  to  the  response  functions  for  normal 
gases  constructed  from  local  conserved  densities.  As  a  consequence,  a  purely  algebraic  manipulation  of  the  response 
functions  leads  to  an  identification  and  representation  of  transport  coefficients  in  terms  of  Green-Kubo  expressions. 

The  Green-Kubo  representation  may  be  a  useful  alternative  form  for  direct  numerical  evaluation  of  the  transport 
coefficients  [16].  More  generally,  the  analysis  given  here  based  on  the  Boltzmann  equation  points  the  direction  for  a 
similar  but  formally  exact  derivation  of  Green-Kubo  relations  from  the  Liouville  equation.  That  analysis  is  in  progress. 
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A.  APPENDIX:  SPECTRUM  OF  RESPONSE  FUNCTIONS 

Using  the  conservation  laws  response  functions  are  found  to  obey  the  equations 

(43) 

(44) 

(45) 


(ds  +  Am(0))G„m(kV)  +ikt  G;im(kV)  =  0 
(dj  +  A„ (0) ) G„m (k* , s)  +  ikjH^m (k* ,  s)  =  8n2  (v,  fc„(kV)) 
(^  +  A„(0))G;im(kV)  +i*J.G&,(kV)  =  8n2  (vX(kV)) 
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where  GJ,m(k*  ,s),  H}lm(k*  ,s),  and  G{‘ln  ( k“ ,  ,v)  are  response  functions  involving  the  fluxes 


<4»(kV)  =  (&» tf„(kV))  =  ($»(kV),)i)  ,  *C.(kV).  =  (%,<MkV))  , 

and  the  source  terms  Snm  (k* ,  s)  and  S'nm  (k* ,  s)  are 


GJ,L(k*,s)=  kV))  , 
(46) 


Snm( k  !■*!  "  ^  (k  !■*))  i  (k  ,  -j)  —  <5„2  (tL'KjiCk  ,■?))  • 


(47) 


Define  the  Laplace  transform  of  a  function  F(s )  by 


F(z)  =  [  dse~zsF(s).  (48) 

Jo 

Then  the  Laplace  transformed  equations  (43)  -  (44)are  (in  an  obvious  matrix  notation  where  I  is  the  identity  matrix 
and  A(0)  is  a  diagonal  matrix  with  elements  A„  ((.))) 

G(k*,z)  (z  +  A(0))  +  ik*G'{k* ,z)  =  /  (49) 

:(z+A(0))G(k*z)  +  iA*^'(k*,z)  =I  +  S(k*,z)  (50) 

(z  +  A(0))G(k*,z)  +  ik*j&i(k*,z)  =  (&,)£,)  +  Si(k*,z)  (51) 

In  this  notation  the  matrix  M(k*,z)  defining  the  dispersion  relations  is  defined  by 

G(k*,z)  =  [(z/  +  A(0))+M(k*,z)]-1. 

Equation  (49)  then  gives 

M(k*,z)  =  /fc,-G-1  (k*  ,z)G'(k*  ,z)  ■  (52) 

Multiplying  by(z  +  A  (0))  and  its  inverse  gives 


M(k*,z) 


ik*  ((z+A(0))G(k*,z))  '  ((z  + A(0)))  G'(k*,z) 

ik*  (f  —  itf H‘ (k* ,z)  +  S(k*,z)^ (-»}G>,'(k*jz)+  ($,?)  +  S'(k*,z)) 


(53) 


where  Eqs.  (50)  and  (51)  have  been  used  in  the  second  equality.  This  is  the  desired  form  expressing  M( k  *,z)  in  terms 
of  flux  response  functions  and  with  the  dominant  k*  dependence  explicit. 

Further  simplifications  occur  for  small  k* .  If  only  terms  through  k*2  are  to  be  retained  (Navier-Stokes  order)  then 
only  //'(O.z)  is  require  in  (53).  Furthermore  5'(k*,z)  is  of  order  k*2  from  symmetry  and  does  not  contribute  to  this 
order,  and  5'(0,z)  =  0.  Therefore 


M(k*,z)  -J.  ik, :  +  ktkj  ^G-,'(0,z)  —  H-i(0,z)  (*,/)  +  (S'(k,z))^0j  .  (54) 

Writing  explicitly  the  transforms  on  the  right  side  as  time  integrals,  and  using  the  identity 

G&.(0„s)  -Hip(0, s)  (jpym)  =  ,  (55) 

gives  Eqs.  (39)  -  (41)  of  the  text. 
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